Renormalization group theory for the phase field crystal equation 
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I. INTRODUCTION 

A fundamental theoretical and computational chal- 
lenge in materials modeling is that of simultaneously cap- 
turing dynamics occurring over a wide range of length 
and time scales, under processing conditions. A clas- 
sic example of such a multiscale problem is dendritic 
growth, a phenomenon seen in the solidification of un- 
dercooled melts, which involves the capillary length ~ 
10~^m, the scale of the pattern ~ 10~^m, and the diffu- 
sion length ~ 10~'*m, which together span length scales 
over five orders of magnitude, and heat/solute trans- 
port through diffusion which occurs on time scales of 
^ 10~^s. Only recently, after considerable advances 
in computing technology, and through the use of so- 
phisticated computational techniques^ 0, 01) has this 
problem become tractable in three dimensions. Al- 
though a number of computational approaches ,4 M 
eluding quasi-continuum methods 0,0,13,01 the hetero- 
geneous multiscale method [TollTll |. multi-scale molecular 
dynamics 12, 13, l3,^3, multi grid variant sflfij and ex- 
tensions of the phase field model T7] have been proposed, 
they all appear to have significant limitations. Develop- 
ing a "handshake" algorithm, to seamlessly integrate the 
transition between scales, is a common problem. Also, 
most models appear to be capable of handling only a 
limited number of crystallographic orientations, with a 
few isolated defects. 

Elder et al. 0, 0| recently proposed a continuum, 
non-linear partial differential equation, which they called 
the phase field crystal (PFC) model, for realistically 
describing materials processing phenomena in polycrys- 
talline materials. The PFC equation describes the evo- 
lution of the time averaged density field in the mate- 
rial, subjected to the essential constraint of mass con- 
servation. While averaging the density makes it possible 
for the model to capture phenomena over diffusive time 



scales (not possible with molecular dynamics), the spa- 
tially periodic variations, with a wavelength on the order 
of the inter-atomic spacing, allow it to incorporate lattice 
defects such as vacancies. These rapid spatial variations, 
however, are also the bane of the model, as they neces- 
sitate the use of a uniformly fine computational mesh 
whose grid spacing approaches the interatomic separa- 
tion. Thus the PFC model is computationally expen- 
sive for mesoscale problems (such as dendritic growth), 
although still considerably better than a molecular dy- 
namics calculation. 

We have recentlydescribed a theoretical approach to 
this difficulty |23, l2ll . presenting a heuristic renormal- 
ization group (RG) method to coarse-grain the PFC 
equation and obtain equations of motion for the ampli- 
tude and phase of the periodic density field. Using these 
variables, it is possible to reconstruct the original field to 
a certain order of approximation. The main advantage of 
such a description is that the amplitude and phase of the 
density field vary on length scales much larger in com- 
parison to those of the field itself, which enables us to 
use a coarser mesh to speed up calculations (see |23| for 
details on accuracy and speedup). Furthermore, the rel- 
ative uniformity of these variables permits solution of the 
equations governing them on an adaptive grid which, we 
anticipate, will further improve computational efhciency. 
We emphasise that the RG procedure is more than a 
naive coarse-graining in real space of the density, because 
it uses the dynamics inherent in the underlying equation 
to project out the long- wavelength, small-frequency be- 
havior. The essential effect of this distinction is that 
valuable aspects of the phase field crystal model, such as 
its native inclusion of elasticity, are preserved. In poll2l| 
we verified explicitly, by numerical calculation, how the 
RG equations are able to accomplish this, and thus are 
suitable for dealing with polycrystalline systems. 

The main purpose of the present article is to present 
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full details of the systematic calculation to derive such 
coarse-grained equations from the PFC equation. A sec- 
ondary goal is to compare and contrast the variety of 
techniques that are available to derive coarse-grained 
equations of motion. A rather surprising finding of our 
work was that when we followed naively the "cookbook 
recipe" for each method, our results were not identical, 
with the RG methods yielding a form of the amplitude 
equation, slightly different from that derived by the clas- 
sical method of multiple scales. The PFC equation obeys 
a local conservation law, and while this by itself can lead 
to a variety of interesting features [2^ , it also brings to the 
fore an ambiguity with the usual implementation of the 
renormalization procedure, something that is not unique 
to conservation laws. This ambiguity is essentially an op- 
erator ordering one, and can be remedied in a straight- 
forward way. Once done, all methods yield the same am- 
plitude equation, even though the technical details are 
quite distinct in the different methods. As a pedagogi- 
cal exercise, we present the analysis for the Van der Pol 
oscillator in the Appendix, once again obtaining consis- 
tent results from all methods when using the approach 
described herein. Our main conclusion is that the renor- 
malization group method is still considerably easier to 
implement than competing approaches, and in particular 
requires no knowledge of the scaling relationship between 
space and time while achieving full rotational covariance 
at lowest order in perturbation theory. 

In the remainder of this Introduction, we review the 
main conceptual developments leading up to the tech- 
niques described in this paper. It is now relatively stan- 
dard in non-linear pattern formation problems to use am- 
plitude equations to uncover universal features of pat- 
tern forming systems. The formalism, first introduced 
by Newell, Whitehead, and Segel (NWS) [H El to de- 
scribe periodic patterns in Rayleigh-Benard convection, 
offers a way to extract the spatio-temporal envelope of 
these patterns, which then allows one to predict the dy- 
namics qualitatively with very little information about 
microscopic details. Unfortunately however, the NWS 
equation, as originally constructed, can only describe the 
dynamics of patterns oriented along the same fixed direc- 
tion, everywhere in space, whereas physical systems often 
produce complex mosaics of patterns with no particular 
orientational preference. Such mosaics arise in real sys- 
tems which are invariant under rotations, and hence any 
equation which is used to study them must also have 
the crucial property of rotational covariance, something 
that is lacking in the NWS equation. Equations with an 
orientational bias can be very difficult to implement nu- 
merically, especially on systems with arbitrarily oriented 
patterns. Nevertheless, the NWS equation embodies the 
important notion of coarse-graining, which has played 
a significant role in shaping the modern day theory of 
pattern formation ^26j . and also forms the basis of our 
approach to multiscale modeling with the PFC equation. 

Gunaratne et al. first derived a rotationally co- 
variant form of the NWS equation using the method 



of multiple scales [2^ |2^ IboI IbH , where they assumed 
isotropic scaling of the spatial variables. They showed 
that the spatial operator in the NWS equation could be 
symmetrized, by systematically extending the calculation 
to higher orders in the perturbation parameter e, the re- 
duced Rayleigh number. They explained that the finite 
truncation of the perturbation series destroyed the rota- 
tional symmetry of the operator, which could however 
be recovered at a higher order. Another important con- 
clusion of their work was that the qualitative behavior 
of pattern formation remained unchanged if one ignored 
higher order corrections, provided the equation itself was 
rotationally covariant. A drawback of their calculation, 
however, was (as with any application of the method of 
multiple scales), the need to guess a priori, the correct 
scaling of space-time variables. In addition, their cal- 
culation required gradual accumulation of operators and 
terms up to ©(e"*), before a rotationally covariant equa- 
tion emerged. 

A more systematic approach emerged shortly after: 
Chen et al. showed how to perform reductive pertur- 
bation theory using RG methods 32], and obtained the 
NWS equation for the Swift-Hohenberg equationjsl by 
renormalizin g th e leading secular divergences at each or- 
der. Graham [SJI subsequently showed that, in fact, this 
method gave the fully rotational covariant equations, if 
all secular terms are renormalized and a careful choice 
of operator splitting is used. Calculations involving the 
RG typically produce elegant and accurate uniformly 
valid approximations for ordinary differential equations 
(ODEs), starting from simple perturbation series where 
no knowledge of the scaling present in the system is 
exercised J2J. For partial differential equations (PDEs), 
the same approach is successful, but generates a tedious 
number of perturbation terms at higher orders. This dif- 
ficulty arises from the need to explicitly construct secular 
solutions of the highest possible order, at every order in e, 
and from a practical standpoint equals (if not outweighs) 
the advantage of requiring no prior insight into the prob- 
lem. For this reason, calculations employing this method 
for PDEs have rarely gone beyond 0{e). The key advan- 
tage of the RG method, however, is that when carefully 
performed, the calculation yields a rotationally covariant 
amplitude equation at a much lower order in e compared 
to the method of multiple scales, as was shown by Gra- 
ham 113. 

Nozaki et al. [sM Is^ , have developed a more abstract 
version of the perturbative RG for weakly non-linear 
PDEs, called the "proto-RG" scheme. They argue that if 
one is willing to sacrifice some of the purely mechanical 
aspects of the conventional RG by taking cognizance of 
the system's properties, such as those exhibited by the 
governing differential equation, one can obtain a rotation- 
ally covariant amplitude equation to 0{e) without having 
to construct any secular solutions. By computing mini- 
mal particular solutions, usually obtained by a straight- 
forward inspection, one can even obtain 0{e^) corrections 
with only a little more algebra. They illustrated the rela- 
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tive simplicity of this method by deriving the rotationally 
covariant form of the NWS equation to 0(e^), as previ- 
ously derived by Gunaratne et al. |23| . Shiwa jl^l further 
demonstrated the efficacy of this scheme by obtaining the 
well-known Cross- Newell phase equation (38i, , which 
describes phase dynamics of patterns generated by the 
Swift-Hohenberg equation. A drawback of this approach 
is in the selection of the so called "proto-RG" operator, 
which turns out to be non-unique in general. Nozaki et 
al. show, however, that the operator is uniquely spec- 
ified, provided we insist on the lowest order differential 
operator possible. 

As the reader will have no doubt realized from this 
synopsis, several methods and variants exist for deriv- 
ing envelope equations from order parameter equations 
(OPEs) that produce predominantly periodic patterns. 
Although one may argue that some of these methods are 
essentially variants of perturbative RG theory for PDEs, 
they are structurally very different. It is thus a very in- 
structive exercise to compare the defining properties of 
each of these methods in the context of a single micro- 
scopic OPE. We present such a detailed study in this 
paper using the PFC equation, which unlike the Swift- 
Hohenberg model, has not been extensively studied. 

The article is organized as follows. To set the context 
for our work, we briefly introduce the PFC model in sec- 
tion m In section IIIII we present a detailed derivation 
of an amplitude equation from the PFC model using a 
heuristic approach. In section HVI we use the proto-RG 
method to derive the amplitude equation more systemat- 
ically. We attempt to verify these calculations indepen- 
dently in section using the method of multiple scales. 
A 1-D derivation via the conventional RG method is pre- 
sented for completeness in section IVTI We find that while 
the proto-RG and RG results are consistent, they do not 
agree with the other calculations, due to an operator or- 
dering ambiguity not previously noticed. We remedy this 
in lVIII and conclude with some remarks in section IVlIII 



II. THE PHASE FIELD CRYSTAL EQUATION 

The phase field crystal model proposes a sixth or- 
der non-linear partial differential equation for describing 
the space-time evolution of the time-averaged, conserved 
density variable i), of a material. As has been shown 
[TsL ITgj . this equation has the potential to accurately 
model a variety of key materials processing phenomena, 
including heterogeneous nucleation and grain growth, liq- 
uid phase epitaxial growth, ductile fracture mechanics, 
dislocation mechanics [40j | , and plasticity. An important 
feature which differentiates it from another popular con- 
tinuum material model, the phase field model 17], is its 
incorporation of elasticity in the free energy functional 
through terms that guarantee gradients in the equilib- 
rium density field, for certain values of the control pa- 
rameters. Crucial to the construction of this free energy, 
is the observation that elasticity is a natural property of 



a system which is characterized by periodic fields. We 
refer the reader to the exhaustive article by Elder and 
Grant ^ detailed description of the model and its 

applications. 

From the point of view of pattern formation theory, 
the PFC equation is the conserved analog of the simplest 
form of the Swift-Hohenberg equation (with only the cu- 
bic non-linearity), and is written as 

= [{(l + v2)2-e}^ + V;3] +C. (1) 

Here, e is the scaled undercooling, a parameter akin to 
the modified Rayleigh number, controlling the stability of 
the uniform phase ^ (liquid) to the appearance of either 
a periodic striped phase or a periodic hexagonal phase 
(both crystalline solids), and C is the conserved Gaussian 
noise which accounts for thermal fluctuations in the sys- 
tem. A phase-diagram illustrating the phase boundaries 
and coexistence curves in e — -i/i space is given in [T9l | . 

For present purposes we will concern ourselves only 
with the uniform to hexagonal phase transition (but are 
not in any way restricted to it), and will disregard C 
in view of its relative unimportance in describing phase 
transition kinetics. In a single mode approximation, a 
hexagonal pattern is described by 

3 

V'(x,i)=^A,(i)e*'^^--fVi + c.c., (2) 

where ki — /co(~i\/3/2 — j/2), k2 = fcoj and ks = 
fco (1-^/3/2 — j/2) are the reciprocal lattice vectors, /cq is 
the wavenumber of the pattern, i and j are unit vectors 
in the x- and y-directions, Aj are the complex amplitude 
functions, and c. c. denotes complex conjugate. We point 
out that a striped pattern, with the stripes parallel to the 
X axis, can be represented by the same equation above 
with ^1 = ^3 = and A2 7^ 0. As shown in Fig. ^ how- 
ever, the PFC equation describes the evolution of several 
such hexagonally patterned crystals of arbitrary orienta- 
tion, that collide to form grain boundaries. While the 
pattern remains periodic within each crystal, there is a 
break in the periodicity across the boundaries. Equation 
^ can be made to describe such a system by allowing 
the Aj to be spatially varying, i. e. ^^(x, t). Our goal 
is then to derive evolution equations for ^^(x, i), which 
along with Eq. (O can be used to reconstruct V'(x, t) 
in a one-mode approximation. Note that Aj(yi,t) now 
contains information about both the envelope function 
(amplitude modulus), as well as the orientation (phase 
angle) of each grain, but varies on a much larger length 
scale. 

In order to proceed with our analysis, which is essen- 
tially perturbative, we identify e as a small parameter, 
in powers of which we shall expand if) about the steady 
state solution. We point out that the PFC model stip- 
ulates no such restrictions on the value of e, other than 
e > 0, whereas it is natural to expect our treatment to 
restrict the validity of the amplitude equations so derived 
to small values of e <C 1. 
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(a) t = 300 





(b) t = 450 
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III. HEURISTIC RG CALCULATION 

We now present a derivation of the amplitude equation 
from Eq. ^ using linear stability analysis and a short- 
cut motivated by experience. An idea along these lines 
was previously implemented by Bestehorn and Haken 
|4l| to derive an OPE (similar to the Swift-Hohcnberg 
equation) for modeling traveling waves and pulses in 
two-dimensional systems, but not for deriving amplitude 
equations. 

We pose Eq. in a more convenient form by scaling 
"0 by and calling this new variable tp. In this manner, 
Eq. becomes 



dtil: = V^(l + V^)^V + - i^)- 



(3) 



Let us now consider the stability of the uniform phase 
solution ^ to the formation of the hexagonal pattern by 
adding to it a small perturbation ip, so that ijj = ijj + ijj . 
Substituting in Eq. (|2Jl and linearizing about ip we obtain 



dt^p = V2 [e(3V? - 1) + (1 + V2)2] ^ 



(4) 



If -0 is a hexagonal instability in the form given by the 
spatially dependent part of Eq. 10), then using Aj{t) = 
Aoj exp{ujjt), where Aoj are complex constants, and sub- 



stituting in Eq. 
relation 



iQ, we obtain the discrete dispersion 



(5) 



after applying orthogonality conditions. Here ujj predicts 
the growth or decay rate of a hexagonal instability in the 
spatially uniform system. Note that for real values of -tp, 
LOj is always real. Thus, a necessary condition for the 
instability to grow, i. e. for ojj to take on positive values, 
is S'f/'^ — 1 < 0, or equivalently 3-0^ — e < in original 
variables. The most dangerous wave-number is the locus 
mi = /co = l. ^ ^ 

We now consider spatial modulations in the amplitude 
about this preferred wave-number, i. e. 



Aj{t) I — * Ai?j(x,<) = Aoj-e" 



,(Q)tg''Q-x 



(6) 



where Q = Qx '\ + Qy j is a perturbation vector, and Ajij 
is the renormalized amplitude, whose implication will be 
clarified in a following paragraph. Consistent with Eq. 
Q, we can now write the exponent controlling growth 
rate along each lattice vector as 




(c) t = 750 



FIG. 1: Heterogeneous nucleation, crystal growth, and for- 
mation of grain boundaries in a 2-D film from three randomly 
oriented seeds, as simulated by the PFC model. The field 
plotted is the density variable i/)(x, f). Note that the pattern 
is periodic inside each grain. 



co.iQ) = |Q+k,f e(l - 302) _ (1 - IQ + k,| 



(7) 



We now replace the Fourier space variables in the above 
equation by their real space counterparts so that. 



LUj = dt , Qx 
thus obtaining. 



IQ- 



^3 I 



1- 



2ikj ■ V 



-idy 



(8) 
(9) 
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Combining Eqs. and (O, the space-time amplitude 
variations along each lattice vector is given by the sixth 
order linear partial differential equation 

We also need nonlinear terms, which play a vital role in 
pattern dynamics near onset of the instability, to com- 
plement the above set of equations. There are a cou- 
ple of different ways to obtain these terms. One can 
directly look for the nonlinear part in the normal form 
equations 1421 for the dynamics of Aj in a hexagonal ba- 
sis 123, 143L l44j . These equations have been widely used 
to study the dynamics and stability of exactly periodic 
rolls and hexagonal patterns originating from the static 
conducting state in Rayleigh-Benard convection. Alter- 
natively, one can derive these terms to a particular order 
in e through a renormalization group (or multiple scales) 
analysis of the governing differential equation. Here, we 
choose the latter approach, starting from Eq. but 
only going far enough in the RG analysis to identify the 
correct form of the terms. 

We start with a perturbation series in e 

V' = V'o + eVi+eV2 + e^^3 + --- (11) 

where ipo is a steady state solution and ^'^(^^o) ^-re the 
higher order corrections. As we are interested in ampli- 
tude variations in the hexagonal pattern, we pick t/jq to 
be the steady hexagonal solution, i. e. Eq. |(2Jl with Aj{t) 
replaced by Aj{t 00). Substituting in Eq. we 
obtain the following equation at 0{e): 

[dt - V2(l + V2)2] = [dt ~ Cr) ^1 = V2(^3 - V-o), 

(12) 

where 

3 

V^i^l-ijo) - (l-3^2)^^^.gik,.x 
-3Ai(|Ai|2 + 2|A2|' + 21^31') e*^-'' 

-3A2 (21^112 + 1^212+21^312)6*"-^ 

-3A3 (2|Ai|2 + 2|A2|2 + |A3|2)e*^-^ 
-6A*2A*i:e'^' '' - 6AtA*V5e*-^ 
-eAlAl^e'^^ '' + other terms + c.c. (13) 

The superscript denotes complex conjugation. To this 
order, the "other terms" are functions of complex expo- 
nentials that do not lie in the null space of the linear dif- 
ferential operator in Eq. (|12|l . i. e. they are non- resonant 
terms. Therefore, they do not contribute to unbounded 
growth in Tpi. The terms listed in Eq. H13I) arc, how- 
ever, resonant with the operator, and their coefficients 
need to be renormalized in order to bound the solution 
obtained by truncating the perturbation series at 0{e). 
The renormalization procedure allows the amplitude Aj , 
previously constant, to now have space-time variations 



that absorb secular divergences. We assert that the non- 
linear terms in the amplitude equation to 0{e) must be 
the renormalized coefficients of the exponential terms in 
resonance with the differential operator. For example, 
the terms complementing the space-time operator along 
basis vector ki must be 

e(l-3^2)^^^ _ SeARi{\ARi\^ + 2\Aii2\^ + 2\Ar3\^) 
- 6eA*j,2A*j,,i,, (14) 

where the Arj are the renormalized amplitude functions 
(no longer constants). Note that these terms are com- 
pletely identical to those pred icted by normal form the- 
ory for a hexagonal basis |23,|4JI- Combining Eqs. HlOf) 
and H14|l we write the amplitude equation as 

dtAi - -{l-Cu,)Cl^Ai-e{l-3i>^)C^,A, 

+e(l - 3V?)Ai - 3eAi(|Ai|2 + 2\A2\^ + 2|A3|2) 
-GeA^Ali;, (15) 

for lattice vector ki , and permutations thereof for k2 and 
k3 , where we have replaced the variables Arj by Aj . 

We observe that the leading term in Eq. (|14|l is consis- 
tent with the right hand side of Eq. H1U|) , thereby provid- 
ing a natural overlapping link about which to match the 
linear stability and perturbation results. In a more ab- 
stract sense, we draw a parallel between this method and 
the technique of matched asymptotic expansions in sin- 
gular perturbation theory, where inner and outer asymp- 
totic solutions are matched over a common region of va- 
lidity in the solution space, to obtain a globally valid 
solution. This completes our derivation of the amplitude 
equation via a heuristic or "quick and dirty" approach. 
For future reference, we will call Eq. H15|l . the QDRG 
(quick and dirty RG) equation, and the method used 
to obtain it as the QDRG (or heuristic) method. As we 
have already demonstrated the remarkable accuracy with 
which the QDRG equation mimics the PFC equation in a 
previous article 20], we will refrain from presenting any 
new evidence to that effect here. 

To summarize the procedure, we first conducted a lin- 
ear stability analysis of the scaled PFC equation about 
the uniform state to obtain a linear differential operator 
controlling the space-time evolution of the complex am- 
plitude Aj of the hexagonal pattern. We superimposed 
on this dispersion relation periodic modulations of the 
amplitude, and from the dispersion relation in terms of 
these latter modulations, identified the gradient terms in 
the amplitude equation. We then carried out the first 
step in a conventional RG analysis to obtain the form of 
the nonlinear terms that should accompany this differ- 
ential operator, and combining the two results, we wrote 
down the amplitude equation for the hexagonal pattern. 
In this respect, our approach lacks the full mathemati- 
cal rigor of a conventional RG reduction or a multiple 
scales derivation, which gives it a somewhat "dirty" ap- 
pearance. However, we made no assumptions about the 
scaling of the space-time variables in the system, nor did 
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we have to construct any secular solutions so far. We 
will comment on extending this method systematically 
to higher orders in e in the following section. 

IV. PROTO-RENORMALIZATION GROUP 
DERIVATION 

With the proto-RG method, our starting point is 
Eq. © with the perturbation series Eq. Hll(l . Thus, to 
0{e) we obtain Eq. (|12|l . whereas to C(e^) we get 

idt~Cv)^2^yHs^o^i-^i)- (16) 

The structure of Eq. (|12|l allows us to infer that its sim- 
plest particular solution will take the form 

3 3 
3 3 
2 2 
2 

+ + c, (17) 

where 

si = -1 — - J 2 ' ^2 = 1 — - J 2 ' S3 = S2 - si 
r3\/3 -1 ^3V3 -1 

tl = —1 1 — ,t2 = 1 1 — 

2 •'2 2 2 

-73 r5 r\/3 r5 

vi = -r\/3-j*2,V2 =r\/3-r2, (18) 

are non-resonant modes generated by the cubic term. 
Note that we have explicitly denoted the space-time de- 
pendence of the secular coefficients Pij (x, t) , which are 
polynomials in x, y and t, whereas by inspection, the 
other coefficients Qij, Rij, Sij, Tij, Uij and Vij can be 
complex constants. Specifically, Pn satisfies 

{dt - Cv) Piie'^' '' = (l-3V;^)Aie*i-^ 

-3Ai + 2|A2p + 2IA3H e*^-" 

-6A*Al^Pe'^'■'' 

^ {dt + {1 ~ C^JClJ Pn = 

+ 2|A2|2 + 2IA3I') - 6A;Ali, 

= Lk.Pu, (19) 



where Li^. is the "proto-RG" operator for lattice vector 
kj. From the above equation it is quite obvious that 
Pii cannot be constant for any non-trivial solutions, and 
likewise for P12 and P13. 

As Pij are secular, we now renormalize |3^ tp about 
arbitrary regularization points X and T, as in the con- 
ventional RG method, to get 

3 

= ^-t-5]A«,(X,T)e*^-- 

3 

+e5](Py(x,t)-Pi,(X,T))e*^- 

+ ... + C.C. (20) 

where A^j is now the renormalized amplitude that ab- 
sorbs secular divergences. Since ip must be independent 
of these regularization points, we have 

L^f^ = 

^ L^;^ARiiX,T) - eL^;^Pni^,T) + e^L^;^P2iO^,T) 
+e'L^fP3i{X,T) + ... (21) 

after applying orthogonality conditions. This is the gen- 
eral form of the proto-RG equation for weakly nonlinear 
oscillators '3^. L^'^ is the proto-RG operator Lki hi 
Eq. (|19|l . with variables x and t replaced by X and T 
respectively. Changing back from (X, T) (x, t) and 
Anj Aj, and using Eqs. p9|) and H21|) . we can write 
the amplitude equation along lattice vector ki to 0{e) 
explicitly as 

dtAi = -(l-£kJ/:ki^i+e(l-37/i2)Ai-3e.4i(|^i|2 
+2\A2f + 2\A3f)^6eA;Al^, (22) 

with appropriate permutations for A2 and ^3. Note 
that in using Eq. ((T^ . we have replaced Aj by their 
renormalized counterparts Anj as is consistent with the 
proto-RG procedure, before reverting to the former nota- 
tion for amplitude. Upon comparing the two amplitude 
equations obtained so far, Eqs. p5|) and (|22|l . we note 
that the QDRG derived equation carries the extra term 
e(l - 3-ip'^)Ck^Ai. Evidently, the QDRG and the proto- 
RG methods produce different amplitude equations when 
applied to the PFC equation, the extent of this difference 
controlled by the parameter e(l — 3'0^). 

As mentioned earlier, the principal advantage of using 
the proto-RG method is the relative ease with which one 
can progress to higher order calculations. Let us now 
extend this calculation to 0{€^). We need ?Ai in order 
to evaluate the right hand side of Eq. (|16|l . which means 
that we additionally need to evaluate Qij, Rij, Sij, Tij, 
Uij and Vij. The constant values of these terms can be 
determined by inspection. For example, by analogy with 
Pii, we see that Qn must satisfy 

{dt ~ Cv) Qiie^'^' '' = -12 {Aj^ + 2AiA;AI) e'^'^'-'' 
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(5t + (4-£2kJ(/:2ki -3)2)Qn = 
-12 {Aj^P + 2AiA^Al 



Thus P21 satisfies 



(23) 



Unlike Eq. p9|) however, we see that Eq. H23|l permits a 
constant solution for Qu, which in turn is determined to 
be 



1 



{Aj^j + 2AiA;A* 



(24) 



Similarly, constant solutions for the other coefficients are 
Q12 = -^{Al{p + 2AIA2AI) 
Qi3 - -^{Al^p + 2AIA*^A,) 



511 = ~^(^AlA3 + 2Ai4,A*2 + AfA* 

512 = -^(^AiAl + 2A;M3 + AfAl 
^ (^A2Al + 2Ali;A3 + AfA*^ 



Sl3 = - 



Til = - 

T12 = - 

Uii = - 

U12 = - 

Vii = - 

V12 = - 



^1^3 

12 
AlAl 

12 
AiAf 

12 

12 

AlAl 
12 

42 A* 
12 



(25) 



We know that the particular solution to Eq. has the 
form 



^^2 



3 

E 



P2j(x, t)e'''^ -'' + 



+ c. c. 



(26) 



where we have shown only the resonant part of the so- 
lution. The terms on the right hand side of Eq. (|16|l . 
resonant with lattice vector ki, evaluate to 

^ = [(i-3v;2-6(i^ip + |A2p + i^3n)(i-/:kjPii 

-3^2 (1 - /IkJ - &A1A2 (1 - ricj A*2 

-QA1A3 (1 - /:kj ^1*3 - 6^A* (1 - /ikj ^1*2 

-6M; (1 - >CkJ ^1*3 - 6AiA; (1 - C^,)Pi2 
-6A1 A* (1 - /Ik J A3 - 6^2^511 - iAlUii 
-QA^AlTii - QA2A\Vii - QA\AlSii - SAfSu 
-3AfSl2 - GAs^pSl^ - 6AlA;Sl3 - SAjT*, 
-6A2A;QI^ - 6A3A*Qt3 - iAfRii 
-6A2A3Qii-6^Atgn]e*i". 



(28) 



The non-constant terms in $ (terms containing Pij) are 
now ignored |35l l3^ while the remaining terms are de- 
termined from their constant solutions, Eqs. H24(l and 
(|25|l . Thus, using Eq. (|21|) we can write the amplitude 
equation along lattice vector ki to ©(e^) as 

dtAi = -{l-£u,)£l^Ai+e{l-34>^)Ai-3eAi{\Ai\^ 
+2\A2\^ + 2\A3\^) - GeA^Ali^ + ll^4,A\A2A3 
+eV^i(2|Aip-K9|A2p + 9|A3p) 
+ lUm2\Ai\^ + \A2\^ + \A3\MlAl 

+ ^e'AlAfAf + 5e^Ai\Ai\' (IA2P + 1^31') 



A2e'Ai\A2\'\A3\' + ^e^Ai\Ai\^ 



3 

64^ 

.^e'Ai\A2\^ + le'Ai\A3\ 



(29) 



with cyclic permutations for lattice vectors k2 and k3. 

We can in principle extend the QDRG method also 
to higher orders by performing the same steps above, 
until the point where we identify the resonant terms on 
the right hand side of Eq. (|16|l . i.e. $. Combining this 
result with Eqs. pU|) and p4|l we can then obtain the 
amplitude equation Ea. (|29|l . but with an extra term e(l — 
3i,^)U,Ai. 

In summary, both the proto-RG and the QDRG can be 
calculated including terms of 0{e^), and the results dif- 
fer by a small but non-zero term. Which, if any, of these 
calculations is correct? And what is the origin of the 
discrepancy between the two methods? Is the QDRG re- 
sult not to be trusted, being derived heuristically? Faced 
with two seemingly incompatible, although very similar 
results, it is natural to attempt an independent test of 
the analysis, which we did using the standard method of 
multiple scales. This calculation is presented below, but 
owing to technical complications arising from the inter- 
ference of modes and the need to go to sixth order of 
perturbation theory, we found it only feasible to perform 
the calculation for the case of one dimension. Neverthe- 
less, we will see that, in fact, the QDRG result, Eq. l(T^ . 
is more correct. The small discrepancy between this re- 
sult and the proto-RG result is finally resolved in Section 

IvITI 



V. MULTIPLE SCALES DERIVATION 

We now re-derive the amplitude equation using the 
traditional method of multiple scales. As the primary 
purpose of this derivation is to verify the previous cal- 
culations via an independent method, we stick to a one 
dimensional analysis here that considerably simplifies the 
(27) algebra. For convenience we use = e, and write Eq. Q 
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in 1-D as 



(30) 



The basic premise of the multiple scales analysis is that 
while the pattern itself varies on the scale of its wave- 
length (27r//co), its amplitude varies on much larger 
length and time scales. It is then appropriate to intro- 
duce slowly varying arguments 



X ^Sx 



T 



(31) 



for the envelope function A{X^T). This scaling was pre- 
viously applied by Gunaratne et al. j23| to the Swift- 
Hohenberg equation with success (based on the form of 
the discrete dispersion relation), and as the PFC equation 
is essentially a conserved analog of the Swift-Hohenberg 
equation we anticipate that the same scaling holds here. 
Derivatives scale as follows 



dx dx+ Sdx 
dl -> dl + 26dxd^ 

whereas the operator 









such that 








01(1 + dl)' 










Adldt + IQdldl (1 






i2dldl + 89ia, (1 


+ 01) 




lW%dl + 29i (1 + 




= 













(32) 



(33) 



(34) 

We now expand -0 in a perturbation series in S to get 

1/; = + S^/Ji + 5^^2 + 6^iJ3 + .... (35) 

Using Eq. H32|) and the above series, the 6 expansion of 
the nonlinear term in Eq. (|30|l can be written as 

+6[dU3i;liJi - V'l) + 2dxdA^i;l - 0o)] 

(37/^0 V'l +37/^^^2-02) 

+29x9. (30^ Vi - 01) + dxiy^o - V'o)] 
H-<53 [9^(02 + 6000102 + 30^03 - 03) 

-f29x9,(30o0? + 30^02 - 02) 
+91(30201- 01 )]+<S4 [92 (30.202 

+30^002 + 60)O0'103 + 30004 - 0'4) 

+2dxd,{iPl + 600011^2 + 30;203 _ 03) 

+9|(30o0? + 3^-^02 -V^2)] 

+C'((55). (36) 



Substituting Eq. (|35() in Eq. H30|l . and using the scaled 
operators in Eqs. I|32I34II . we can write equations satisfied 
by the 0^, at each 0{S"'). At 0(1) we obtain, 

£o0o = 

=>0o = + Aoi(X,T)e" + c.c. (37) 

where Amn is the complex amplitude of mode n at 0{5™). 
At 0{6) we get 



£o0'i + ^1^0 = 



c.c. 



(38) 



where (and henceforth) wc neglect the constant term in 
view of its inclusion in Eq. 137|l . At the next order we 
have 

£o02 = 9t0o - Ci^i - /:20o - dli^^ - Vo). (39) 

For 02 (a:^, t) to remain bounded we have to guarantee that 
the right hand side of Eq. does not have a projection 
in the null space of Cq, which yields a solvability condi- 
tion j2fij iilj (also known as the Fredholm alternative). 
Applying the alternative imposes the following condition 
on the amplitude at 0{S'^): 

OtAoi = 49^^01 + (1 - 3i'^)Aoi ~ 3^oi|^oi|'. (40) 



Thus, 



1p2^Ml^ +^226 +^236 +C.C 



(41) 



where A22 = ^§^0/3, and A2 



i22 = ^01 V^/-^' anu /i23 = ^oi/64. 
At subsequent orders, the following equations are ob- 
tained for 0m: 



0{5^) : CqiPs, = dripi - Ciip2 - - ^300 

-[92(30201 - 0i) + 29x9,(03 - 0o)] 

OiS"^) : £o'04 = 9t02 - >Ci03 - ^202 - ^^301 - AV'o 
-[92(30^00? + 3V'o'02-V^2) 
+29x9,(31^2^1 - 0i) + d],i^Pl - 0o)] 

0{S^) : £o^J5 = Ot^IJs ~ -^104 - £20'3 - -^302 - AV'l 
-rsV'O - [9'(0'? + 600011^2 + 30;2^3 

-V>3) + 29x9,(3V'o^? + 30^1^2 - 02) 

+9i(3020;i_V;i)] 

-£50)1 - Letpo - [92(3021/^2 + 3i/'o0| 

+ 60)o0'l0'3 + 30'oV'4 - 04) + 29x9,(02 

+60;o0'i0'2 + 3V'o0'3 - V'a) + 5i(30oV^? 
+30;o2^2-V^2)], (42) 



and successive applications of the Fredholm alternative 
yield the following amplitude equations at those respec- 
tive orders. 



OtAu = -12^91^01 +49i^ii 



(1 - 3^^)2idxAoi 



+ (1 - 3i/;2)Aii - 6^ii|Aoi|2 - 3^2^AJi 
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+&idx{Al^Al{) 

dTA2i - -139^Aoi- 12^91^11+49^^21 
-(l-3^2)(2*9xAii+aiAoi) 
+ (1 - 3^2)^21 - SAl^A*^ - 6|AoipA2i 
-6^PA22A*^ ~ 3^23^*1^ - 6Aoi|Aii|2 
-3Al,A*2, + 6tdx{Al,Al, + 2|Aoi|'Aii) 
+SdliAlM 

drAsi = 6z9iAoi - 139^^11 - 12za|-A2i + 4^1^31 
-{I - 3^P^){2idxA2i + dj,An) 

+ (1 - 3^2)^3^ _ 6AiiA2iA*i ~ 6IA01PA31 

-GAoiAaiAJi - 6^23^5171^1 - 6AoiAiiA*i 
-3Al,Al, - 6i'A32A*^ - 3A33A51' 
-6M22AI1 + 6i9x(2Aoi|Aii|2 + 2IA01PA21 

+^01^11 + ^01^21) + (^01^11 
+2|^oi|2A*i)+h.o.t. 

5t^4i = 9^.101 + 6iai All - 139^^21 - 12z9|A3i 

+4aiA4i - (1 - 3,p^){2tdxA3, + aiA2i) 

+ (1 - 3^2)^41 - SAl^A*^ - 6A11A31AS1 
-6Aoiyl4iyl*i - 6^11^21^*1 - 6^01^31^*^ 
-3Al,A*2, - 6Aoi|A2i|2 - 6AoiAiiA*i 
-SAl^Al^ - 6Mi2A*^ - 3A43AS1' 
-6^A32Al^ ~ 6A33A*^Al^ - 3A23AJi' 
-6^A33A;2 - 6A01IA23P + 6idx{2AiiA2iA*^ 
+2IA01PA31 + Aiil^iip + 2AoiA2iAl^ 
+2AoiAiiA;^ + Al^Ay + 3d%{2Aoi\Aii\^ 
+A*^Al^ + 2IA01PA21 + AI^A;^) + h.o.t. (43) 

Here, "h.o.t." refers to higher order terms that arc func- 
tions of Aqi and its derivatives. The amphtude function 
for the pattern (e*^) can be written as 

A{X, T) = Am {X, T) + Mil (X, T) + 5^A2i (X, T) + . . . . 

(44) 

Using Eqs. (gUJ), and (jHll . and scahng back to orig- 
inal variables, i.e. X — > 5~^x and T — > the ampli- 
tude equation to 0{5'^) can be written as 

d^A = AdlA-l2idlA-lidtA + QidlA + dlA 

-6^{1 - 3^2) {2id^ + dl)A + [(1 - 3V7)A 

-3A\A\' + 3 {2^^, + dl) (AlA^)] ^*(^A|A|4 
+24,^A\A\^) + 0{S'') (45) 

or more compactly, after replacing 5^ e, to 0{e^) 

dtA = -il-C,D)C,D''A-eil-'3^^)CiDA 

+e{l - 3iIj^)A - 3eA\A\^ + 3eCiDiA\A\^) 

-e^{^A\Af + 2^'A\A\')+0{e'), (46) 
where Cid = {2idx + d^). 



Let us now compare Eq. H46|l with the one dimensional 
equivalents of the 0{e) amplitude equations that we have 
previously derived for hexagonal patterns, i.e. Eqs. H15|l 
and (|22|) . Without re-deriving, the one dimensional 
equivalents are readily obtained by setting 

^2 = A3 = and £k. = CiD (47) 

in those equations. We observe that Eq. (|46|) (also trun- 
cated to 0{e)) contains at least one term which is not 
present in either of the equations previously derived. 

We note that the QDRG result of Eq. ((TB|) is closer 
to the multiple scales result compared to the proto-RG 
result (and the RG result in sectionlvT)!. in that it fails to 
capture only the nonlinear derivative term at 0(e), which 
is actually a higher order correction to 3yl|ylp. This is 
clearly because the spatial operator in the QDRG method 
is an outcome of a linear stability analysis, whereas one 
would have to perform a nonlinear stability analysis to 
obtain nonlinear spatial derivative terms. The clear ad- 
vantage of the QDRG calculation however is that it was 
done with significantly less effort, and in a rotationally- 
covariant manner; perturbation theory to 0(e) was all 
that was required. The multiple scales analysis, on the 
other hand, required a sixth order perturbation theory 
treatment, and in order to simplify the algebra, we only 
worked in one dimension. In higher dimensions, the in- 
terference between the modes would have created a huge 
increase in the complexity at each successively higher or- 
der in perturbation theory. The QDRG calculation is 
only heuristic, but as we will show below, can be justi- 
fied from a full calculation, albeit with a minor technical 
modification of the previously-published recipe, to take 
into account the special feature of the conservation law 
in the PEG model. 

We conclude therefore that although the QDRG re- 
sult of Eq. pSIl and the multiple scales method to 0{e) 
still do not yield consistent results, the QDRG method is 
still an improvement over the proto-RG method. In or- 
der to track down the source of the discrepancy, we next 
attempted a full RG calculation without any shortcuts, 
i.e. systematically calculating explicitly and renormaliz- 
ing all the divergent terms to 0(e). 

VI. RENORMALIZATION GROUP 
DERIVATION 

In this section, we present a derivation of the amplitude 
equation using the conventional RG method, in one di- 
mension for pedagogical simplicity (just as done for the 
method of multiple scales). The calculation is compli- 
cated because of the need to obtain explicit formulae for 
the secular divergences [3^ . but this is possible at the 
order to which we worked. 

Starting from Eq. (|30|1 (with replaced by e) and a 
naive perturbation series in e as in Eq. pil) the zeroth 
and first order solutions can be written as 

i/jg = ■(/; Ae" + c.c. 
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Pi(x,t)e'^ + Qie^" + i?ie^" +C.C. (48) 



The difficulty in the conventional RG method comes from 
the need to explicitly determine the form of the secular 
coefficient Pi. While this is a routine task for ODEs, it is 
far from trivial for PDEs. A further complication is that 
the solution for Pi must be the highest order polynomial 
that satisfies the PDE, to be able to eliminate all secular 
divergences. It turns out that this is critical to obtaining 
the rotationally covariant operator at a lower order in e. 
Using the method of undetermined coefficients wc find 
such a solution to be 

6 

-e(l-3V?-3|^nA^C,Pi,(a:,i), (49) 



where 



Pii 

Pl2 



= t 



Pl3 - 

Pl5 = 
Pl6 = 



(-89280^2 + 7680i^ + A320ixt - M560ixt^ 

720 

-4680xh + 2880a;^t^ - lUOixh + UOxH + x^) 

— (-5760ii^ - 1560xt + 960xt'^ - 120ix^t 
720^ 

+80x^t + x^) 

1 
312 



(192r - 2SSixt + ASxr + x*) 



--(24x< + x3) 



(50) 



and the constants Cj satisfy X]j=i Q = 1- 

The RG method proceeds as follows: (1) dummy vari- 
ables X and T are introduced, (2) the divergent terms in 
Pij of the form a;"*" are split to read = (a;™^ - 

X™T") + X™T", (3) the constant ampHtude A is rede- 
fined using an e expansion A ~ Aii[X, P)(l+X]j=i ^''^j)^ 
where An is now the renormalized amplitude, and Zj are 
the renormalization constants which are chosen order by 
order in e to absorb the X™r" terms, and (4) since the 
solution -0 is independent of X and T, all derivatives of 
^p with respect X, T, or a combination thereof must be 
zero. This last condition yields the following RG equa- 
tions at 0(e), 

dAR 
dT 

dX^ 
O^Ar 
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I2i 



dX^ 
d^A 



R 



dX^ 
S^A 



R 



dX^ 



= Cie(l 


-37/^2 


-?,\A\^)A 


= C2e(l 


-3i/^2 


-3|A|2)A 


= C3e(l 


-3i/^2 


-3|A|2)A 


= C4e(l 


-37/^2 


-?,\A\^)A 


= C5e(l 


-3i/i2 


-3|A|2)A 


= C6e(l 




-i\A\^)A. 



Further, using X]j=i — ^ ^"^^ replacing Ar — > A, 
X ^ X, and T t, the above equations can be combined 
to read 

dtA + (1 - Cid)Cid^A = e(l - 3^^) A - ieA\A\^ (52) 

which is also the 1-D proto-RG equation. 

We close this section with some interesting observa- 
tions, (i) The equations in H51(l do not form a unique set 
of solvability conditions. Other equations are possible, 
e.g. 



1 O^Ar 
16 dX^dT^ 
ii d^AR 
""SdXdT^ 



C2e{l-3ij^ -3\A\^)A 
C3e(l-3V?-3|A|2)A 



(53) 



(51) 



The choice of Eqn. H51|) is motivated by the observation 
that it yields a rotationally covariant amplitude equation, 
and other physical considerations such as the microscopic 
equation being only first order in time, (ii) The list of 
possible terms Pi^ does not include the leading polyno- 
mial term Bx, where B is an arbitrary constant, as this 
term is annihilated by the kernel of the PDE. Thus no 
constraint is available to fix P. It turns out that unless 
this term is also renormalized, all secular divergences are 
not removed. This may explain the absence of certain 
terms in Eq. H52(l that however show up in the multiple 
scales analysis. To be certain, the calculation needs to 
be carried out to higher orders; but we do not attempt 
this here. 



VII. OPERATOR ORDERING AMBIGUITY 
AND ITS RESOLUTION IN THE RG METHOD 

In this section, we resolve the discrepancy between the 
answers generated by the QDRG method, the RG meth- 
ods, and the method of multiple scales. Curiously, no 
such discrepancy was observed in the treatment of the 
Swift-Hohenberg equation, a non-conservative OPE, by 
RG methods [33, yj, |3g and multiple scales techniques 
[27l |. In fact, it can also be readily ascertained that the 
QDRG method will produce the same result as the other 
methods for this equation, which we leave as a simple ex- 
ercise for the reader. Why then does a discrepancy arise 
in the PEG equation? Glearly, the role played by the ex- 
tra Laplacian, a consequence of the conservation law in 
this case, must be non-trivial! 

Note that this Laplacian operator carries over to the 
right hand side of both Eq. H12(l . the 0{e) equation for 
the RG methods, and Eq. (|39|l . the 0{e) equation for 
multiple scales. However also note that, in the method of 
multiple scales, in addition to the non-linear terms, this 
operator is also subjected to an e expansion. There is no 
provision in any of the RG methods to allow the same to 
happen to the Laplacian. In other words, the operator 
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may very well have not existed on the right hand side 
at 0{e), and we would have obtained exactly the same 
result as before! 

A clue to the subtlety is to look at the way in which 
the secular terms are renormalized. The naive way, as 
followed here, would be to evaluate the right hand side 
first, look for secular terms later, and then renormalize 
these divergent coefficients. However, this will not elimi- 
nate secular terms generated by the differential operator. 
In order to eliminate all secular terms, the amplitude 
must be renormalized before differentiation, for the sim- 
ple reason that renormalization and differentiation are 
non-commutable operations. In other words, there is an 
operator ordering ambiguity in the implementation of the 
renormalization group method, exposed in this problem 
by the conservation law. Performing the calculation with 
the operations of renormalization and differentiation re- 
versed is equivalent to performing an e expansion in the 
differential operator. 

We find that by following this procedure, additional 
terms in the coefficients of the resonant modes are auto- 
matically generated. Specifically, when we evaluate the 
right hand side of Eq. (|12|) after assuming the amplitudes 
of "00 to have a space-time dependence, the renormalized 
coefficients of the resonant exp(iki • x) forcing term work 
out to be 

e[(l - 3V?Mi - 3Ai + + 

-6A;Ali;-{l-S^i}^)Ck,Ai 

+6 (l^ip + + £kiAi + 6AJ|VAi|2 

+3AlCi,^Al + 6AiA*Ck,A2 + ^A^A^L^^A^ 

+%AyAt,L^,Ai + %A^A^L^,,Al -f G^^rki^; 

-h6V;A;£ki^3 + 12^Vyl^ • VA^ 

4^12^1 (V^i • V^t + VAs • V^2 + VA3 • \IAl) 

+12A2VA1 ■ VA*2 + 12A*2VAi ■ VA2 

+i2A3VAi-va; + i2a;vAi-va3], (54) 

which when specialized for the 1-D case becomes 
e[(l - 3i/^2)^ _ 3^|^|2 _ (1 _ s^-^)CidA 



+6\A\^CidA + 'SA^CidA* + 6A* 
dA dA\ 



OA 
dx 



+12A ^ ^ , 
ox ox 

= e[(l - 3ij'^)A - 3A\A\^ - (1 - SV^^id^ 

+3Cid{A\A\^)]. (55) 

We note that the above terms are identical to the 0{e) 
terms on the right hand side of Eq. (|46|l . Therefore the 
correct amplitude equation to 0{e) should contain all 
the terms in Eq. 1)54(1 . In order to illustrate the gener- 
ality of this approach, we apply this idea again in the 
appendix to the Van der Pol oscillator, another equation 
for which the previously reported implementation of the 
RG method, and the method of multiple scales produce 
different answers. 



We wish to point out that the assumption of a con- 
stant amplitude in the ipQ solution makes it possible for 
the coefficients of the non-resonant terms in tpi to assume 
constant values, a fact that is favorably used in extend- 
ing the proto-RG calculation to the next order. However, 
with our modification to the proto-RG procedure, it is 
clear that for the PFC equation at least, non-resonant 
coefficients cannot have constant values. Thus comput- 
ing higher order corrections to the amplitude equations, 
will require explicit construction of particular solutions, 
which may limit progress beyond 0{e) by purely analyt- 
ical methods. 



VIII. CONCLUSION 

We have presented a detailed illustration of various 
perturbative techniques to derive amplitude equations 
from order parameter equations that produce periodic 
patterns. Amplitude equations serve as powerful analyt- 
ical tools with which to investigate pattern stability and 
defect interactions, as well as accurate coarse-grained de- 
scriptions of pattern forming systems, and this calls for 
practical and reliable mathematical methods for deriving 
them. 

Although our benchmark for accuracy is the widely- 
accepted method of multiple scales, it is critical to note 
that this method is not failsafe, because it requires a pri- 
ori identification of the way in which space and time scale 
with the small parameter e. There are many instances 
where surprising scales emerge that would not easily be 
identified a priori (e.g. see the analysis of the Mathieu 
equation in [321). 

The method of multiple scales typically involves a very 
lengthy calculation before a rotationally covariant oper- 
ator ensues, and involves computation of various higher 
order terms which ultimately do not improve the over- 
all result significantly. In the example presented here, a 
sixth order calculation was required to get the lowest or- 
der amplitude equation. The reader should bear in mind 
that the fairly involved calculation shown in this article 
was only one dimensional. 

On the other hand, the practicality of RG based meth- 
ods, where the amplitude equation was obtained very 
quickly at 0{e) itself, is self-evident. No guesswork was 
required to determine the scaling of the variables and 
all calculations started with naive perturbation expan- 
sions in e. In particular, our so called "quick and dirty" 
(QDRG) method and the proto-RG method are attrac- 
tive techniques, because there is virtually no need to 
construct explicit solutions. Both methods use only in- 
formation available from the differential equation and in 
that sense, are very general ways of building a controlled 
coarse-grained approximation to the order parameter 
equation being studied. Furthermore, the QDRG method 
gives the correct result quickly, apart from a small non- 
linear rotationally covariant gradient term which is not 
captured by the linear stability argument. 
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At 0{e), we have shown that the QDRG method pro- 
duces a more accurate amphtude equation compared to 
the proto-RG method, by capturing certain extra terms 
that are revealed in the multiple scales analysis. How- 
ever, with our corrected order of operators in the way in 
which the RG is implemented, we find that all methods 
converge identically at this order. 

In conclusion, we have presented a detailed calculation 
of the coarse-graining of the phase field crystal equation, 
for small e. Elsewhere |23,|23|, we have demonstrated the 
utility of the coarse-grained equation in performing large- 
scale simulations of materials processing phenomena in 
two dimensions. Further developments of these tech- 
niques are underway, implementing adaptive mesh refine- 
ment to solve the amplitude equations obtained here in 
two and three dimensions, and going beyond the single 
mode approximation which has formed the basis of our 
work and that of Elder and collaborators. We hope to 
report on these developments at a future date. 
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APPENDIX: VAN DER POL OSCILLATOR 

In this appendix, we explore the commuting of dif- 
ferentiation and renormalization with a simple ordinary 
differential equation example: the Van der Pol oscillator. 
Note that this yet another case in which a differential op- 
erator is multiplied by a small parameter (see right hand 
side of Eq.^lJ). 

The autonomous ODE is given by 



we obtain 



y" + y = eil-y^)y' 



(A.l) 



where ' denotes differentiation with respect to the vari- 
able t. As there is a derivative on the right hand side of 
this equation we anticipate that the proto-RG amplitude 
equation will fail to capture certain terms that turn out 
in the multiple scales analysis. 

It is known that the scaling t — et works for this prob- 
lem Isfl. Hence, 



y' {dt + edr)y 

y" ^ {dj + 2edrdt + ^dl)y, 



(A.2) 



where the subscripts denote partial differentiation. Ex- 
panding 2/ in a perturbation series 



0{1) 
Oie) 



id? + l)yo 
id! + l)yi 

i)y2 







id! 



-2drdtyo + (1 - yDdtyo 
-2drdtyi + (1 - yo)dtyi - d'^yo 
-2yoyidtyo + (1 - yo)dryo- (A.4) 



From this we find 

yo 
yi 



Aai{T)e'* + c.c. 
An(T)e^*+Ai3(T)e3'*-)-c.c 



(A.5) 



Application of the Fredholm alternative at 0(e) and 
0{e^) yields the following amplitude equations 



01 



2idrAoi = (1- I^OlP) 

2idrAii = i(An-2An|^oi|'-^oi^n) 



-a, (Aoi-Aoil^oiP) 
^01 1 A 



on 



(A.6) 



which can be combined after scaling back to original vari- 
ables to get 

d^A + 2idtA = e [iAil ~ \Af) + dtA{l - \A\^)] + ©(e^). 

(A.7) 

Nozaki and Oono 36] on the other hand have obtained 
the following equation using the proto-RG method 



d^A + 2idtA = eiA{l ~ \A\^) + 0{e^). 



(A.8) 



Note that the missing term dtA{l — | Ap) can be captured 
by differentiating the lower order result, i. e. 



2idtA = eiA{l - {Al"^) 



(A.9) 



but this does not seem a very general approach. In par- 
ticular, it is not obvious how this can be extended to 
PDEs. 

The Oie) equation using the proto-RG method reads 



where 



2/1+2/1 = (l-?/o)2/o> 



yQ — Ae'* + c.C. 

yi = P(/)e" + Qe=^** + c.C. 



(A.IO) 



(A.ll) 



where A can be a constant while P cannot. Thus, the 
proto-RG operator turns out to be 



C = d! + 2idt, 
and the proto-RG equation reads 

CA^eCP + 0{e^), 



(A.12) 



(A.13) 



y = 2/0 + £2/1 + e 2/2 + 



(A.3) 



where A is now the renormalized amplitude. When eval- 
uating CP however, we allow for the possibility that A, 
which appears on the right hand side of the equation can 
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also be a function of t, or equivalently renormalize A on 
the right hand side before differentiating j/qi which gives 
us 

LP - e [iA{l - \Af ) + dtA{l - \Af)] . (A.14) 
Therefore the true amphtude equation should read (using 



Eq.HHl) 

£A = e [iA{l - + dtA{l ~ \A\^)] + 0{e^), (A.15) 

which is identical to the multiple scales result of 
Eq. iXTIl . 
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